(N 



> 
O 



Cauchy problem with data on a characteristic 
cone for the Einstein- Vlasov equations* 



^ ; Yvonne Choquet-Bruhat Piotr T. Chrusciel 
§ : March 6, 2013 

u 

O" 1 Dedicated to Mario Novello 



1 Introduction. 

In recent papers (henceforth denoted by I and II) we have proved exis- 
tence and uniqueness theorems for solutions of the Cauchy problem for the 

O 



Einstein equations in vacuum with data on a characteristic cone Co, with 
vertex at O, generalising a construction of [5l[7]. We have used the tensorial 
splitting of the Ricci tensor of a Lorentzian metric gona manifold V as the 
sum of a quasidiagonal hyperbolic system acting on g and a linear first order 
operator acting on a vector H, called the wave-gauge vector. The vector H 
vanishes if g is in wave gauge; that is, if the identity map is a wave map from 
(V,g) onto (V,g), with g some given metric, which we have chosen to be 
Minkowski. The data needed for the reduced PDEs is the trace g of g on Co, 
but the geometric initial data is the degenerate quadratic form g induced by 
g on Co, the missing part of g is determined by a hierarchical system of or- 
dinary differential equation^], called the wave-map-gauge constraints, along 
the rays of Co, deduced from the contraction of the Einstein tensor with a 
tangent to the rays (it is also possible to prescribe g up to a conformal factor 
and impose a coordinate condition). 



*Vienna preprint UWThPh-2012-21 

1 For previous writing of these equations in the case of two intersecting surfaces in 
four-dimensional spacetime see Rendall [7] and Damour-Schmidt [S|. 
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The above generalises in a straightforward manner to Einstein equations 
coupled to matter fields whose energy-momentum tensor satisfies specific 
structure conditions. This is the case for Maxwell fields and scalar fields. 
However, it is not immediately clear that kinematic matter sources described 
by a Vlasov distribution field / fit into this scheme. In fact, there is a basic 
issue arising, related to the fact that the support of the distribution function 
/ has to be contained within the subset of the tangent bundle where the 
momentum of the particles is timelike future pointing; outside of this region 
the Vlasov particles are tachyons. But a no-tachyons condition requires 
a priori knowledge of the whole metric on the initial data hypersurface, 
while in the wave-map scheme above only part of the metric is known before 
the constraints are solved, the remaining components being determined by 
ODEs along the generators of the hypersurface. It has been shown in [I] 
how to modify the scheme so that the whole metric can be prescribed on 
the initial characteristic hypersurface, which solves that problem, and allows 
e.g. the treatment of Vlasov particles with a given mass. However, the 
formulation in I and II has the clear advantage, that only geometrically 
significant objects are prescribed as free data. This is not the case with [I], as 
many components of the initial metric have a gauge character. It is therefore 
of interest to see whether kinetic matter fits into the original scheme of [7], 
as generalised in I and II. In this article we show that this is indeed the case. 
More precisely, we prove that the Cauchy problem on a characteristic cone 
for the Einstein- Vlasov system can also be split into a hierarchical system 
of ordinary differential equations as constraints and an evolution problem, in 
the important case of astrophysical studies where the masses of "particles" 
take a continuous set of positive values. We show how to construct physically 
relevant initial values in this case. 



2 Einstein- Vlasov system. 
2.1 Distribution function. 

In kinetic theory the matter is composed of a collection of particles whose 
size is negligible at the considered scale: rarefied gases in the laboratory, 
galaxies or even clusters of galaxies at the cosmo logical scale. The number 
of particles is so great and their motion so chaotic that it is impossible to 
observe their individual motions. It is assumed that the state of the matter 
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in a spacetime (V, g) is represented^] by a "one particle distribution function": 

Definition 1 A distribution function is a scalar function f > on the 
tangent bundle TV to the spacetime V, 

f : TV —> R by (m,p) M- f(m,p m ), with m£V, p G T m V. 

The physical meaning of the distribution function is that it gives a mea 
number of particles with momentum p G T m V at a point m EV. 

In astrophysics the particles are stars or even galaxies, they do not take 
a finite collection of a priori given masses, but rather a continuous family of 
positive masses. 

One denotes by uj g and u p respectively the volume forms on V and T m V: 

9 = Ulg A Up 

where, in local coordinates 

uj g = (det g)^dx° A dx 1 A ... A dx n , u p := (det g)^dp° A dp 1 ... A dp n . 

Definition 2 The moment of order k of the distribution function f is the 
symmetric k-tensor 

T ai - a "(m) := [ f(m,p)p ai ...p a ?Lu p . 

JT m V 

The moment of order zero, integral on the fiber T m , m G V of the distri- 
bution function: 

M(m) := / fu p 

JT m V 

gives the density of presence of particles at a point m G V . The second 
moment of / defines the stress energy tensor it generates. It is a symmetric 
2-tensor on V obtained at any point m G V by integrating on T m the product 
of f(m,p) by the tensor product p ®p. 

2 A mathematical justification of the oncoming of chaos in relativistic dynamics is a 
largely open problem. 

3 In the sense of Gibbs ensemble. 
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2.2 Vlasov equation. 



In a Lorentzian spacetime (V,g), in the absence of non gravitational forces 
and collisions, each particle follows a geodesic of the spacetime metric g, i.e. 
an orbit of the vector field X = (p, P) on TV, 

P°:=^, |£ = ^> with P a :=-r>V 

with A a canonical parameter and T"^ the Christoffel symbols of the metric g. 
Recall that the scalar g(p,p) is constant along an orbit of X in TV, recall also 
that (Liouville theorem) the volume form 9 is invariant under the geodesic 
flow; that is, Cx denoting the Lie derivative with respect to X, it holds that 

C x 9 = 0. 

In a collisionless model the physical law of conservation of particles im- 
poses to the form f9 to be invariant under the vector field X, hence f9 has a 
zero Lie derivative with respect to X. Since we already know that Cx9 = 
the invariance reduces to the Vlasov equation for / 

df df 

which says that in the phase space TV the derivative of the distribution 
function in the direction of X vanishes 
A fundamental theorem is 



Theorem 3 If the distribution function f satisfies the Vlasov equation, its 
moments have vanishing divergence. In particular 

V„T Q/3 = 0. 



2.3 The Einstein- Vlasov system. 

Given a smooth manifold V, the unknowns are the metric g on V and the 
distribution function / on TV. They must satisfy the Einstein- Vlasov system 

S a/ 3 = T a p, on V, C x f = on TV, 
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where S a p is the Einstein tensor of g, hence satisfies the contracted Bianchi 
identities, V a S a/3 = and 

T a p := / PaP/3f(-,p)uJp 

is the second moment of the distribution function /, hence such that when 
the Vlasov equation holds 

V a T a/3 = 0. 

The Einstein- Vlasov system is therefore coherent. 

3 Cauchy problem on a characteristic cone 

The proof of local existence for solutions of the Cauchy problem with initial 
data for g on a spacelike manifold S and for / on TS can be found in [1] 
and references therein. In what follows we consider the case where the initial 
manifold is a characteristic cone of g. 

The future characteristic cone Co of vertex O for a Lorentzian metric 
g is the set covered by future directed null geodesies issued from O. We 
choose as in I and II coordinates y a such that the coordinates of O are 
y a = and the components g Xfl (0,0) take the diagonal Minkowskian values, 
(— 1, 1, . . . , 1). If ^ is Lorentzian and C 2 in a neighbourhood U of O there is 
an eventually smaller neighbourhood of O, still denoted U, such that Co H U 
is an n dimensional manifold, different iable except at O, and there exist in 
U coordinates y := (y a ) = (y°, y l , % — 1, . . . , n) in which Co is represented 
by the equation of a Minkowskian cone with vertex O, 

C : = {r - y° = 0}, r := 

The null rays of Co are represented by the generators of the Minkowskian 
cone with tangent vector I. We overline as in I and II traces on Co, and 
underline components in the y coordinates. The components of t are (r_ = 1, 
P_ = r~ l y l . We denote by Y Q the causal future of O, y° > r. We denote by 
TCo the tangent bundle with base Co, by T^Co the subbundle with fiber 
future timelike vectors for the metric g. 
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3.1 Cauchy problem for the Vlasov equation with data 
on Co, given g. 

Let the n + 1 dimensional spacetime (V, g) be a given Lorentzian manifold. 
The Vlasov equation, with X = X(g), and unknown /, 

C x f = 0, 

is a first order linear partial differential equation for the distribution function 
/ on the tangent bundle TV, differential equation on the geodesic flow which 
reads in the coordinates y, p of TV 

df(y(X),p(X)) 
dX 

where y and p are solutions of the differential system 
dy a dp a 

~dx ~ v - ' d\~- 

This quasilinear first order differential system has, if its coefficients are Lip- 
shitzian i.e. if the considered metric g is C 1,1 , one and only one solution for 
A — A small enough, (y l (X), y°(X),p a (X)), taking for A = A given values 

y%X ) = y\ y°(X ) = r, p a (A ) = 7r Q , 

7r a , a vector in R n+1 timelike for g. The image in TVof the obtained curves 

X^(y\X), y\X),f(X)), 

is the set of future timelike geodesies (we mean the curves and their tangents) 
issued from points of T + Co in a neighbourhood of O. Well known proper- 
ties of timelike geodesies (in particular the conservation of g(p,p) under the 
geodesic flow) show that conversely a past directed timelike geodesic issued 
from a point (m,p) G T + Yq, with m in a small enough neighbourhood of O, 
will meet T + Co- 

The initial data for / is a function / on T + Co- The solution of the 
Cauchy problem for the Vlasov equation with this initial data is given in a 
neighborhood of T + in T + Yo by 

f(y(X),p(X)) = f(y a (X ),p a (X )). 
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3.2 Cauchy problem for g with data on Co for the Ein- 
stein equations in wave gauge, given the stress en- 
ergy tensor of /. 

The Einstein tensor in wave gauge forms a quasilinear system of wave oper- 
ators deduced from the identities 

<](<?) = -\g X »dl^ + q^(g)(dg,dg). 

The stress energy tensor of a given distribution function fx on a spacetime 
(V, <?i) is an integral operator on f\ and g\ 

Ti >a p = I Pl,aPl,pfl(;P)^ gi - 

One considers the system of quasilinear equations for a metric g2, with T\ 
known, 

Einstein w (£ 2 ) = H. 

If T\ is smooth there exists (Cagnac - Dossa theorem [6]) to > such that 
these equations have a solution g 2 with trace on Co H {r G [0, to]} the trace 
g of a smooth spacetime metric. 

3.3 Solution of the Einstein-Vlasov reduced system 

One shows by iteration using known theorems that the Einstein-Vlasov sys- 
tem in wave gauge with unknowns g and / admits a solution with initial 
data g and / if g is the trace of a sufficiently smooth spacetime tensor field 
and / is the trace of a sufficiently smooth function on TV, compactly sup- 
ported on each future timelike cone, with for example / = and g trace of 
the Minkowski metric in a neighbourhood of the vertex. 

3.4 Solution of the original system 

To show that a solution of the reduced system satisfies the original equations, 
if the initial data satisfy appropriate constraints, we proceed as in I where 
we have proved the decomposition £ a S a/ 3 = C a + C a , where C a depends 
only on g and C a is a linear homogeneous differential operator on the wave 
gauge vector. The Bianchi identities show, as in I, that a solution of the 
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reduced system satisfies the original equations if the initial data satisfy the 
constraints, namely 

We now show that in our astrophysical setting these constraints form again 
a hierarchical system of ordinary differential equations if we use coordinates 
x a adapted to the null structure of Co, defined by 

n n 1 , A A , 1 A\ 

x = r — y , x = r and x — /i (r y ), 

A = 2, ...n, local coordinates on the sphere S n_1 , or angular polar pseudo 
coordinates, such that the curves x° = 0, x A =constant are the null geodesies 
issued from O with tangent £ := -J^. The trace g on Co of the, a priori 
general, Lorentzian metric g that we are going to construct is such that 
g u = and giA = 0; we use the notation 

g = g 00 (dx ) 2 + 2u dx°dx l + 2u A dx°dx A + g AB dx A dx B . 

In particular for the Minkowski metric 77 = 77 it holds that 

^00 = -!> u o = V A = 0, f] AB = r 2 s AB , 

where s AB dx A dx B is the metric of the round sphere S 71 ^ 1 . The volume ele- 
ment in the tangent space to V at a point of Co is in the ^-coordinates 

u! p = vq det g. 

4 Geometric initial data and gauge functions. 

The geometric data is a degenerate quadratic form g induced on Co by our 
unknown g\ it reads in coordinates x 1 , x A 

g = g AB dx A dx B , 

i.e. gn = giA = while g AB dx A dx B = g AB dx A dx B is an x 1 -dependent Rie- 
mannian metric on S' n_1 . While g is intrinsically defined, it is not so for ^ 0) 
isq, is A , they are gauge-dependent quantities, the unknowns of our constraints. 

As in the spacelike case, the second fundamental form x of Co appears in 
the constraints, but here it depends only on the induced metric because the 
normal to Co, £ == a^r, is also tangent. We have defined x i n 1 as the tensor 
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X on Co given by the Lie derivative C^g with respect to £ of the degenerate 
quadratic form g, namely in the coordinates x 1 ,x A : 

Xab = 2 9 i9ab, Xai = Xn = 0. 
We denote by r the mean extrinsic curvature r := g AB XAB- 

5 The first constraint operator 

We have found in I the identity 

FS W = R n ee -dxr + v^vor - ^(T, + r) - x A B X £ A 
By definition of the wave-gauge vector H we have 

I\ ee W 1 + Hx, W 1 = -v,x l g AB s AB , 
and hence Rn decomposes as 

Rn ee Ci + L\ , where £1 := —-Hit vanishes in wave gauge, 

while C\ involves only the values of the coefficients g AB and vq on the light- 
cone; the first constraint reads 

Ci := -Sir + {^o - \{Wx + r)}r - X a B X B A = f u = (^o) 2 T 00 , 

f 00 = detg [ (p°) 2 f(.,p)d n+1 p, d n+1 p = dp°dp 1 ...dp n . 

The first constraint contains uq as only unknown if g and / are given. The 
resulting ODE is singular for r = since r multiplies OiUq. However the al- 
ternative used by Rendall [7] and Damour and Schmidt [5] of prescribing only 
the conformal class of g, splitting x i n t° its conformally invariant traceless 
part a and its unknown trace r and imposing to z/ to annul the parenthesis 
in C\ does not work here because T n depends upon z/ . 

We look for a solution of the first constraint tending to 1 as r tends to 
zero, hence set Uo = 1 + u and write the first constraint as the first order 
differential equation (recall that u° = (iso)" 1 ) 

rdiU = r{a + b) + f3 + (a + 2b + 3/3)n + (b + 3(3)u 2 + (3u 3 , 
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where 



To avoid analytical difficulties near the vertex, we assume that there exists a 
neighborhood < r < r of the vertex O within Co on which the initial data 
/ vanishes and g coincides with fj. In this neighbourhood the coefficients 
reduce to their Minkowskian values 

n — 1 n — 1 

«o = —z — ; »o := ~ — i Po = 0> 

It 2r 

the solution u tending to zero as r tends to zero is then u = for r < r . 
For r > r it is the solution w vanishing for r = r of a smooth equation with 
smooth coefficients, equivalently of the integral equation 

u(r,x A )= [ {(a + b) + (x l y 1 [/3 + (a + 2b + 3p)u + (b + 3P)u 2 + fiu 3 ]}dx l . 



This solution can be computed by the Picard iteration method as long as 
the interval [ro, r] is small enough; its size, and the norm of u depend on the 
size of the various coefficients. We will always have z/o > in a small enough 
neighbourhood, as required for the Lorentzian character of g. 

6 The Ca constraint 

We have written in I the Ca constraint operator. The Ca constraint with 
kinetic source reads 

C A - T\a = - l -m A + ri A ) + V bxI ~ \d A r + <9a(^ + u^) - f 1A = 0, 
where £ A is defined as 

U := -lifidvA + Au% cX C a + (w° - v A + 9ab9 CD {S B d - t B CD ), 



while 

Tia(x) 



[ p lP AF(x,p)d n+1 p = is { [ P°(vaP° + g A BV B )d n+1 p. 



When z/ has been determined the only unknown in Ca~T\a is va- It satisfies 
a linear equation, its solution follows the same lines as in I. 
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7 The Co constraint 

When vq and v A have been determined, the Cq constraint operator, £PSop = 
Soi, contains as only unknown g 00 = 9oo, we have found in I and II after 
long computations that the Cq constraint can be written, using the other 
constraints 

5iC+(^ + r)C+^{9 1 W 1 + ( K + r)W 1 + i?-^^ + ^V^B} + ^/ = 0, 

h f :=^g n T u + g 1A f 1A + g 10 f w , 

k = v Q d x v Q --(Wx + t), Wi = is W°, W° = W 1 = -rg AB s AB . 
The relation between ( and 500 is 

where g 11 and (foo are linked by the linear relation 

9oo = 9oo = -9 {v Q ) + g v B v A . 

The component T 01 of the stress energy tensor of the distribution function 
reads 

foi(x) = / PoPif(x,p)d n+1 p = u (g 00 p° + v Q p l + v A p A )p°J{x,p) d n+l p. 

The only unknown remaining in the Cq constraint is goo- It appears linearly. 
The solution equal to —1 in [0, r ] is smooth and negative in [0, r + e\. 

8 The tachyon problem. 

For a physical reason, namely the positivity of the masses of the consid- 
ered "particles" , the initial data / of the distribution function must have its 
support in the subset 

g o{P°T + 2zW + 2u A p°p A + g AB p A p B < 0. (1) 

Since only g AB has been a priori given the above inequality gives an a posteri- 
ori restriction on the support of /. It is natural to assume that this support in 



11 



momentum space at the tip, where all components of the metric are known, 
is contained within the future Minkowski timelike-cone. Then, by continu- 
ity of the data / and of the solutions of the constraints, one obtains some 
neighborhood of the tip in Co where this support will remain in the future 
cone of the constructed metric. This neighborhood will provide the relevant 
initial data. A precise statement to this effect can be obtained as follows: 

Recall that p° = p° + p 1 is the time component of p in the y coordinates, 
the condition that we will give is that p is strictly timelike for a known 
metric on the support of /. The following lemma is a result of the continuity 
properties found for g: 

Lemma 4 Assume that there exists r > such that the support of the data 
f in TCo is contained in the subset 

(p ) 2 > k{g AB p A p B + (p 1 ) 2 }, k > 1 and p° > (2) 

with g~Ab — Vab f or r — r o (possibly r = 0). Then there is a neighbourhood 
r < r < R in Co such that the support of f in TCo H { r o < r < R} includes 
no tachyons for the metric g. 

Proof. The no-tachyon (positive masses) condition can be written with 
vq = 1 + u, g 00 = — 1 + a where u, a and v& tend to zero as r tends to r 
(possibly zero) 

(p0)2 _ (pl)2 _ a(p )2 + 2p0{up l + ^ _ - ABp A p B > Q (3) 

Remark that p° + p 1 = p° implies 

(p ) 2 < 2{(p ) 2 + (p 1 ) 2 }, 2pV < (P°) 2 + (P 1 ) 2 < 2(p ) 2 + Sip 1 ) 2 , (4) 
while setting \v\ := (J) Ab vaVb)^ we have 

WaP A \ < (9abP A P B )^\v\ hence 2p°v A p A < \v\{(p ) 2 + g AB p A p B }- (5) 

The inequalities (4) and (5) show that the no-tachyon condition (3) is implied 
by the following inequality 

(p°) 2 -(2|«|+4| W | + |^|)(p°) 2 > (pl)2 + ^ p A^ +3 |^| (p l ) 2 + |~| {(p ) 2 + - Ai5p A p i 3 

(6) 
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The assumption (2) of the lemma (which depends only on the given data 
g) implies that there exist two numbers k± < 1 and k 2 > 1 with k^ l k 2 = k 
such that 

k 1 (p°) 2 >k 2 {g AB p A p B + (p 1 ) 2 }, k>l andp°>0. (7) 

The properties of a, it and v, which tend to zero as r tends to ro, show that 
there exists R > r such that for r < r < R it holds that 

(p0)2 _ (2 | a | + 4 | u | + |^|)(p0)2 > fci( p 0)2j 

(p 1 ) 2 + g AB p A p B + 3Mb 1 ) 2 + \v\g A BP A p B < k 2 {g AB p A p B + (p 1 ) 2 }. 
These inequalities complete the proof of the lemma. ■ 



9 Conclusion 

From a smooth solution of the constraints in the x variables, Minkowskian 
for < r < r , one deduces that g is the trace on Co of a smooth spacetime 
function, which completes the proof of the existence of a solution of an 
Einstein- Vlasov spacetime in a neighbourhood of O taking the given data 
g, f. Geometric uniqueness is also easy to prove. The existence up to the 
vertex is more involved: the technique of admissible series used in II could 
be used, but requires further work. 
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